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Abstract 


This  note  is  a  continuation  of  previous  work  on  the  singularity  expansion 
method.  Different  integral-equation  formulations  describing  electromagnetic 

t  l 

scattering  from  imperfectly  conducting  bodies  are  considered.  A  set  of  volume- 
surface  integral  equations  is  used  to  determine  the  analytical  properties  in 
the  complex  frequency  plane  of  the  field  scattered  from  imperfectly  conducting, 
finite  bodies.  Conditions  are  determined  for  the  constitutive  parameters  of 
the  scattering  body  so  that  the  scattered  field  can  be  described  only  by  damped 
sinusoidal  oscillations  when  the  incident  field  is  a  delta  function  plane  wave. 

Scattering  from  a  perfectly  conducting,  finite  body  within  a  parallel 
plate  region  is  also  considered.  It  is  shown  that  the  singularities  in  the 
complex  frequency  plane  of  the  scattered  field  are  poles  and  branch  cuts.  The 
locations  of  the  branch  cuts  depend  only  on  the  separation  between  the  parallel 
plates. 
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I.  Introduction 

The  singularity  expansion  method  for  solving  electromagnetic  interaction 
problems  is  first  introduced  in  Interaction  Note  88^^.  Scattering  from  a 
perfectly  conducting  body  of  finite  extent  is  considered  in  Interaction  Note 
92l  .  It  is  shown  in  C23  that  the  operator  inverse  to  the  integral  operator 

of  the  magnetic  field  formulation  is  an  analytic  operator-valued  function  in 
the  complex  frequency  (s)  plane  except  at  certain  points  in  the  left  half  plane 
where  it  has  poles.  A  representation  of  the  inverse  operator  in  terms  of  the 
natural  modes  and  frequencies  is  also  given  in  reference  [2]].  A  scheme  for 
the  numerical  evaluation  of  the  natural  modes  and  frequencies  for  perfectly 
conducting  bodies  of  revolution  has  been  developed  in  reference  C273*  In  E273 
the  magnetic  field  integral  equation  is  simplified  to  account  for  rotational 
symmetry.  A  computer  code  for  the  numerical  evaluation  of  the  natural  frequencies 
and  modes  is  also  presented  in  C 2713  together  with  detailed  numerical  calculations 
of  some  of  the  natural  frequencies  and  modes  for  a  perfectly  conducting  prolate 
spheroid . 

The  natural  frequencies  and  modes  of  a  thin  cylinder  has  been  calculated 

numerically  in  reference  [253  by  using  the  electric  field  Integral  equation. 

Some  approximate  analytical  results  of  the  locations  of  some  of  the  natural 

frequencies  of  a  thin  cylinder  have  been  derived  in  reference  [263.  The 

integral  equation  is  solved  approximately  by  using  Fourier  transform  methods 

combined  with  the  Wiener-Hopf  technique.  The  solution  thus  obtained  gives  an 

accurate  prediction  of  the  locations  of  the  natural  frequencies  of  perfectly 

("253 

conducting  cylinders  with  diameter-to-length  ratio  less  than  1/100  .  However, 

due  to  the  approximations  introduced,  the  method  fails  to  give  the  correct 
analytical  properties  in  the  complex  frequency  plane  of  the  field  scattered 
from  a  thin  wire.  For  example,  the  approximations  introduce  additional  branch 
points  which  should  not  be  present. 

In  this  note  we  will  consider  the  general  problems  of  scattering  from 
imperfectly  conducting,  inhomogeneous  bodies  in  free  space  and  of  scattering 
from  perfectly  conducting  bodies  within  a  parallel  plate  region.  The  mathematical 
methods  that  we  use  resemble  in  a  way  those  used  in  reference  [23. 
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In  section  III,  electromagnetic  scattering  from  an  imperfectly  conducting 
body  is  first  formulated  in  terms  of  a  volume  integral  equation  for  the  fields 
inside  the  body.  The  mathematical  formalism  used  in  section  III  is  similar  to 

[9-121 

the  one  used  in  deriving  the  so-called  Oseen  extinction  theorem  .  The 

integral  equation  is  of  the  second  kind.  However,  due  to  the  singularity  in 
the  kernel  of  the  integral  equation  we  are  unable  to  draw  any  conclusions  about 
the  analytical  properties  of  the  solution  of  the  integral  equation.  Therefore, 
in  section  IV,  we  reformulate  the  electromagnetic  scattering  problem  in  terms 
of  a  set  of  volume-surface  integral  equations.  When  the  scattering  body  is  a 
pure  dielectric  body  this  set  of  integral  equations  reduces  to  an  integral 
equation  previously  used  in  connection  with  considerations  of  uniqueness  and 
existence  of  the  solution  of  certain  electromagnetic  scattering  problems^*^. 

This  set  of  volume-surface  integral  equations  is  solved  by  using  the 
Fredholm  determinant  theory.  From  this  solution  it  follows  that  the  scattered 
field  has  two  types  of  singularities  in  the  complex  frequency  plane.  The  first 
type  is  due  to  the  singularities  of  the  incident  field  (waveform  singularities). 

The  second  type  is  due  to  the  scattering  body  itself.  Sufficient  conditions 
for  this  latter  type  of  singularities  to  consist  of  only  poles  are  (1)  the 
scattering  body  is  of  finite  extent  and  (2)  the  constitutive  parameters  o,  t, 
p  together  with  Ve'/e'  and  Vp/p  are  analytic  functions  in  the  entire  complex 
frequency  plane.  Here,  e'  =  £  +  io/w,  and  a,  e,  p  are  the  conductivity, 
permittivity  and  permeability,  respectively,  of  the  scattering  body. 

A  knowledge  of  the  order  of  each  pole  is  necessary  when  solving  transient 
electromagnetic  interaction  problems  by  using  the  singularity  expansion  method. 

As  of  now,  no  general  analytical  method  of  determining  the  order  of  each  pole 
has  been  found.  Therefore,  the  order  of  each  pole  has  to  be  determined  numerically 
when  solving  a  specific  problem.  However,  a  general  method  for  the  numerical 
evaluation  of  the  order  of  each  pole  has  been  discussed  in  reference  C2 7 J . 

Scattering  from  a  perfectly  conducting  cylinder  within  a  parallel  plate 
region  has  been  treated  in  references  [203  through  [24]  by  using  the  conventional 
method  of  first  solving  the  integral  equation  numerically  in  the  frequency  domain 
and  then  applying  a  numerical  inverse  Fourier  transform.  In  section  V  of  this 
note  we  will  investigate  the  analytical  properties  of  the  field  scattered  from 
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a  perfectly  conducting,  finite  body,  located  within  a  parallel  plate  region. 

The  effect  of  the  parallel  plates  can  be  taken  into  account  by  using  the  method 
of  images.  Due  to  the  interaction  between  the  body  and  the  two  parallel  plates 
we  will  get  an  infinite  set  of  image  bodies.  Therefore,  the  theory  developed 
in  reference  C2]  cannot  be  applied  directly. 

The  series  representation  of  the  Green's  function,  derived  from  the  method 
of  images,  converges  only  in  the  right  half  of  the  complex  frequency  plane.  In 
order  to  make  use  of  the  singularity  expansion  method  we  first  have  to  analytically 
continue  the  Green's  function  into  the  entire  complex  frequency  plane.  This  is 
done  by  finding  an  integral  representation  of  the  series  defining  the  Green's 
function.  From  this  integral  representation  it  follows  that  the  kernel  of  the 
integral  equation  has  branch  cuts  in  the  left  half  plane.  The  locations  of  these 
branch  cuts  are  uniquely  determined  by  the  distance  separating  the  parallel  plates. 

Once  a  representation  of  the  Green's  function,  valid  in  the  entire  complex 
frequency  plane,  is  available,  it  is  easy  to  show  that  the  operator  inverse  to 
the  integral  operator  of  the  magnetic  field  formulation  has  two  types  of  singular¬ 
ities.  One  type  is  poles  at  those  frequencies  where  there  exist  nontrivial 
solutions  of  the  homogeneous  integral  equation.  The  other  type  is  branch  cuts 
coinciding  with  the  branch  cuts  of  the  Green's  function. 

Before  starting  with  the  analysis  outlined  here  we  will  in  section  II 
briefly  consider  a  surface  integral  equation  describing  scattering  from  a 
homogeneous,  imperfectly  conducting  body. 
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II.  Scattering  From  a  Finite,  Imperfectly  Conducting,  Homogeneous  Body 

Before  we  consider  scattering  from  an  inhomogeneous,  imperfectly  conducting 
body  we  will  discuss,  in  this  section,  briefly  the  special  case  of  scattering 
from  a  homogeneous,  imperfectly  conducting  body.  Mathematically  this  scattering 
problem  can  be  formulated  in  terms  of  a  set  of  two  coupled  surface  integral 
equations. 

Let  the  scattering  body  occupy  the  region  V  of  finite  extent.  The  boundary 
surface  of  V  is  denoted  by  S  and  the  outward  unit  normal  of  S  is  denoted  by  n. 

We  also  assume  that  the  electromagnetic  properties  of  the  scattering  body  can 
be  described  by  the  constitutive  parameters  o,  e,  y.  These  parameters  are  allowed 
to  vary  with  frequency,  but  they  do  not  vary  with  position  within  V.  (See 
figure  1).  Throughout  this  section  we  use  harmonic  time  dependence  exp(-iwt). 

Introduce  the  effective  electric  surface  current  and  the  effective 
magnetic  surface  current  k,  defined  by 


i  =  H*H 

(2.1) 

k  *  nxE. 

(2.2) 

Assuming  that  all  sources  of  the  electromagnetic  field  are  outside  V  we  can 


derive  the  following  set  of  two  coupled  integral  equations  for  J[  and  k 


DD 


\  i  -  k-i  -  z;1  M-k  -  iinc 
1 JE. +  V*  -  zc^*i  -  0 


(2.3) 


where 


n(Ox[VG(£,r’  ;ik)x£(r.’)IldS' » 

JS 


(2.4) 
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(M-£)(r)  =  (ik)"1  j  n(r)x{k2  G(r,r ' ;ik)£(r ' ) 

3S 

~  C<u(r')  +  Kv(r')lTC(r,r'  ;ik)In(r')  -  n(r)D*£(r) 

-  [&(!')  -  g(r)]-^7G(r,r';ik)|dS',  (2.5) 

where 

G(£,£*;ik)  =  (4tt  |r_  -  r_'  |  )-1  exp(ik|ir  -  jr '  | ) 

the  free-space  Green's  function,  V'  is  the  tangential  gradient  with  respect  to  r', 

s  f  41  f 

and  <v  are  the  principal  curvatures  of  S  ,  and  /  indicates  the  principal 
value  of  the  integral,  i.e.,  ^ 

]  • * •  =  lim  • • • 

JS  6-*0  JS-SX 

O 

Here  S-S^  is  the  part  of  the  surface  S  outside  a  sphere  with  radius  6  and  center 

at  _r.  The  subscript  c  in  (2.3)  on  the  operators  L^,  and  indicates  that  the 

wave  number  k  of  the  medium  should  be  used  in  the  Green's  function  instead 
c 

of  the  free-space  wave  number  k.  Also, 

k£  =  wCu(t  +  io/oOJ4,  (2.6) 

Z£  =  Qi/(e  +  io/u))^4,  (2.7) 

Z£  being  the  wave  impedance  of  the  medium,  Zq  «  377G,  and  c  is  the  vacuum  speed 
of  light. 

We  will  now  investigate  some  of  the  properties  of  the  integral  operators 
and  ]£  defined  by  equations  (2.4)  and  (2.5),  respectively.  In  reference  [2] 
it  has  been  shown  that  the  kernel  of  the  integral  expression  (2.4)  behaves  like 

i  t  ~1  2 

|jf  “  i  I  as  jr'  ■+•  _r.  From  this  fact  it  follows  that  the  operator  jy.  ,  defined 
by  ■  Jk’(lt*i),  is  °f  Hilbert-Schmidt  type^"2^. 

In  order  to  investigate  the  properties  of  the  operator  £1  we  assume  that 
<u  and  <v  exist  everywhere  on  S.  Let  P  be  a  point  on  S  with  coordinates 
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r  *=  (x,y,z).  (See  figure  2).  We  choose  the  coordinate  system  in  such  a  way 
that,  in  the  vicinity  of  P,  the  surface  S  can  be  described  in  the  following 
way: 

S  -  |  r '  *  (x'.y'.z')  :  x'  -  x  *  p  cos  ip, 

y*  -  y  ■  P  Sin  ♦,**-*■  p2f(<|0  +  0(p3)}.  (2.8) 

Here  0  £  ip  £  2ir  and  0  <  p  £  e,  where  ei<u  <<  1  and  e*v  <<  1.  The  representation 
(2.8)  is  valid  everywhere  on  S  since  S  is  a  smooth  surface.  Assuming  that  the 
function  £(r)  is  differentiable  on  S  we  then  have 

£(r')  -  £(r)  -  P£(1)(r)cos  ip  +  p£^(r)sin  ip  +  0(p2)  (2.9) 

where 

£(1)(r)  -  (9/3x)*(r), 

£(2)(r)  -  (3/3y)£(r). 

Moreover,  for  an  arbitrary  vector  a  we  have 

(a.V)VG  -  af(R)  +  R(a-R)[R_1  f'(R)3  (2.10) 


where 


f  (R)  -  (4irR)-3(ikR  -  l)exp(ikR),  R  -  r  -  r’ 


and  the  prime  denotes  differentiation  with  respect  to  R.  After  some  algebraic 
manipulations  we  get 


[&(r'>  -  s/r^'^VGCr.jr' ; ik)  ■  (4*p2)  (r)cos  *l>  +  g^(r)sin  ip 

-  3|g^(r)cos2ifi  +  gy2^  (r)sin2i|> 

+  cos  sin  i|£gyl)(r)  +  g^(r)Dj 

Cx  cos  ip  +  y  sin  t{i3  |  +  0(p  *).  (2.11) 


In  order  to  evaluate  the  integral  (2.5)  we  proceed  as  follows.  Let 
be  the  part  of  S  inside  a  sphere  with  radius  e  and  center  at  it.  Moreover,  let 
us  denote  the  integrand  of  the  integral  expression  (2.5)  by  H(r,r')  so  that 

(M-&)(r)  *  f  H(r,r')dS'.  (2.12) 

”  J  S 


We  then  have 


Since  e 
when  jr ' 
def ined 


f  H(r,r’ 
'  S— S 


)dS* . 


(2.13) 


(M-^)(r)  -  f  H(r,r')dS'  + 

JS 

e  t 

is  finite  it  is  easy  to  see  from  equation  (2.5)  that  )  is  finite 

belongs  to  S-SE-  Thus,  the  last  integral  of  equation  (2.13)  is  a  well- 
function  of  r.  From  equation  (2.11)  it  follows  that 


H(r,r’)dS'  *  lim  f  (4ir)-1[N(iJ»)  -£(1)  (r)  +  £<♦>  -&(2)  (r)]p-2  dS’ 

6+0  >S  -S, 

€  e  6 


Here  (£,£*) 


+  lim 
6+0  J 


S 


-S 

e 


j^Cr.r'MS'. 

6 


(2.14) 


as  r 


£.  Therefore,  the  principal  value  Integral, 


lim  f  H. (r,r')dS', 

6+0  'S  -Sx 

can  be  replaced  by  an  ordinary  integral  over  S^.  Furthermore,  £(t|>)  and  £(i|») 
are  two  matrices  with  elements 
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n  (tp) 

*  cos  ip  —  3 

3, 

cos  ip 

XX 

n  (ip) 

2 

*  -  3  cos  ip 

sin  ip 

xy 

n  (ip) 

*  n  (4>) 

yx 

xy 

(2.15) 


n  (ip)  *  cos  ip  -  3  cos  \p  sin 

yy 


q  (<p)  ■*  sin  -  3  cos  ip  sin  ip 

XX 


q  (^)  =  _  3  cos  ip  sin  ip 
xy 


q  (ip)  »  q  (ip) 
Myx  Hxy 


q^(ip)  «  sin  ip  -  3  sin  ip. 


(2.16) 


From  the  expressions  (2.15)  and  (2.16)  it  follows  that 

e  2it 


t  r  2x 

N(tp)dip  *  £0p)dip  *  0. 
*  n 


(2.17) 


Thus,  we  have 


lira 


(4tt)~1CN=(tP)  *^(l>  <tL)  +  g(<P)-&(2)(r)]p ~2  d S'  =  0.  (2.18) 


6-+0  JS  -S. 
e  6 

Equation  (2.18)  implies  that  the  operator  is  a  bounded  operator  when  operating 
on  differentiable  functions  and  the  surface  of  integration  is  a  smooth  surface. 
From  equations  (2.14)  through  (2.16)  it  follows  that 


|H(r,r')|  ~  | r  -  r '  | 


.1-2 


as  r  -*■  r 


(2.19) 


and  that  the  integral 
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I  |H(r,r') jdS'  (2.20) 

H 

e 

does  not  exist.  From  equation  (2.19)  it  follows  that  the  kernel  of  any  iterated 
operator  $J.n,  defined  by 


Mn.£  =  M*  (Mn_1  .£> ,  n  >  1 
_2 

also  behaves  like  |r_  -  r_'  |  as  _r'  "*■  (c.f.  reference  C 19 D> .  Thus,  there 

exists  no  n  such  that  the  operator  £tn  is  of  Hilbert-Schmidt  type.  Therefore, 
the  Fredholm  determinant  theory  cannot  be  used  when  solving  the  set  of  integral 
equations  (2.3). 

Next,  we  will  study  the  analytical  behavior,  in  the  complex  frequency  plane, 
of  the  operators  L  ,  M ,  M  .  By  making  the  substitution  y  =  -  ik,  k  ■  iy 
where  y  is  any  complex  number  we  can  extend  the  validity  of  the  set  of  integral 
equations  (2.3)  to  all  complex  frequencies  s  =  cy.  From  equations  (2.4)  and  (2.5) 
it  easily  follows  that  L  and  |i  are  analytic  operator-valued  functions  of  y  in 
the  entire  complex  frequency  plane.  It  also  follows  from  equations  (2.3)  through 
(2.5)  that  the  integral  operator  defined  by  the  left  hand  side  of  the  set  of 
integral  equations  (2.3)  is  an  analytic  operator-valued  function  of  y  provided 
that  Yc  and  Z c  are  analytic  functions  of  y.  Here 

Yc  =  c{yey[y  +  o/Cce)]}*5  (2.21) 

Zc  =  ^MG  1  y^y  +  a/(cE)3}li.  (2.22) 

Unfortunately,  the  fact  that  o,  t  and  y  are  all  analytic  functions  does  not 
necessarily  ensure  that  yc  and  Z are  analytic  functions  of  y.  For  example, 
when  o,  e  and  p  are  constants  we  notice  from  equations  (2.21)  and  (2.22)  that 
Yc  and  Zc  have  branch  Points  at  y1  =  0  and  Y2  “  “  o/(ce).  Thus,  the  operators 
and  ^  are  analytic  operators  in  the  entire  y-plane  except  along  the  branch 
line  C,  defined  by 


C  =  {y  :  -  o/(ce)  s  Re{y }  S  0,  Im{y)  =  0} . 
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In  sections  III  and  IV  we  will  study  electromagnetic  scattering  from  an 
inhomogeneous, imperfectly  conducting  body.  When  o,  e,  u  and  the  incident  field 
are  analytic  functions  of  y  we  will  show  that  the  only  singularities  of  the 
induced  volume  current  density  are  poles.  This  fact  implies  that  the  solution 
of  the  set  of  integral  equations  (2.3)  is  an  analytic  function  of  y  except 
for  poles  when  a,  e,  y  and  the  incident  field  are  analytic  functions  of  y. 
Unfortunately  we  can  not  draw  this  conclusion  from  the  set  of  integral  equations 

(2.3)  for  two  reasons.  The  first  reason  stems  from  the  fact  that  Z  and  y 

c  c 

are  not  necessarily  analytic  functions  of  y  even  if  c,  e  and  y  are  analytic 
functions  of  y.  The  second  reason  is  due  to  the  fact  that  the  Fredholm 
determinant  theory  can  not  be  used  when  solving  the  set  of  integral  equations 

(2.3) .  However,  in  view  of  the  results  in  section  IV  the  set  of  integral 
equations  (2.3)  can  be  used  in  finding  the  natural  frequencies  and  modes  of 

a  homogeneous,  imperfectly  conducting,  finite  body.  It  is  a  well-known  fact 
that  a  surface  integral  equation  is  better  suited  for  numerical  treatment 
than  a  volume  integral  equation. 

In  the  next  section  we  will  derive  a  volume  integral  equation  for 
calculating  the  scattering  from  an  inhomogeneous  imperfectly  conducting  body 
of  finite  size. 
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III.  Volume  Integral  Equation  Describing  Scattering  From  an  f 

Inhomogeneous,  Imperfectly  Conducting  Body 

In  this  section  we  will  derive  a  volume  integral  equation  that  describes 
electromagnetic  scattering  from  a  body  whose  conductivity  (a) ,  permittivity 
(e),  and  permeability  (u)  are  functions  of  space  and  frequency.  Throughout 
this  section  we  will  suppress  a  harmonic  time  dependence  factor  exp(-iwt). 

From  Maxwell's  equations  one  can  derive  the  following  set  of  partial 
differential  equations  for  the  electromagnetic  potentials  ($,4)  in  an  arbitrary 
medium 


V24>  +  k2$>  =  -  e_1(p  -  V.p) 
o  — 

(3.1) 

V2A  +  k2A  =  -  UQ(i  +  VxM  -  iwP) . 

Here  k  =  w/c,  P  is  the  polarization  of  the  medium,  M  is  the  magnetization  of 
the  medium,  and  p  and  i^  represent  the  charge  and  current  density,  respectively, 
of  exterior  sources.  The  electromagnetic  fields  £,  ji  are  obtained  from 


E  =  iu)A  -  V$ 


B  =  V*A. 


(3.2) 


We  also  have  the  relationships 


E  -  e"1 (D  -  P) 

B  -  y  (H  +  M) 
—  o  —  — 


(3.3) 


where 

V-D  -  p 

VxH  -  _i  -  iwD. 

From  equations  (3.1)  and  (3.4)  it  follows  that 


(3.4) 
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V-A 


ike  =  0. 


(3.5) 


Next,  we  consider  a  region  V  of  finite  extent  with  some  polarization  P 
and  magnetization  M.  Futhermore,  we  assume  that  there  are  no  exterior  sources 
of  the  electromagnetic  field  in  V  implying  that  p  =  0  and  i^  =  0  in  V.  Let  S 
be  the  boundary  surface  of  V  and  let  n  be  the  outward  unit  normal  of  S.  The 
complement  of  V  with  respect  to  the  three  dimensional  Euclidean  space  is  denoted 
by  Vj.  We  also  assume  that  the  region  is  vacuum  as  far  as  the  electromagnetic 
properties  are  concerned.  Assuming  that  all  the  exterior  sources  of  the  electro¬ 
magnetic  field  are  at  infinity  we  have  the  following  differential  equations  for 
the  electromagnetic  potentials 


2  2-1 
V^<t>  +  k  *  =  e  v*P 
o  — 

V2A  +  k2A  =  yQ(iwP  -  VxM)  in  V 


(3.6) 


and 

2  2 
V  +  k  <t>  *  0 

V2A  +  k2A  =0  in  Vj . 


(3.7) 


Let  G  be  the  free-space  Green's  function,  and  let  r  be  a  point  inside  V. 

We  then  have 

(GV2$  -  <J>V '  2G)dV '  -  V *  •  (GVt>  -  4>V'G)dV ' 

Jv-v6  Jv-v6 

r.(GV4>  -  4>V'G)dS'  (3.8) 

where  V'  operates  on  the  second  argument  of  G(r,r_' ;ik),  n'  *  n(r'),  and  is  a 
sphere  with  surface  S^,  radius  S,  and  center  at  r_  (see  figure  3).  It  is  easy  to 
show  that 


n'.(GV4>  -  $V'G)dS 


» 

I  _ 

. 


lim  R’  (GV4>  -  4>V'G)dS'  -  4>.  (3.9) 

6-*0  Js6 
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Equations  (3.6),  (3.8)  and  (3.9)  then  give 


4>  +  e"1  GV.PdV'  +  (4>3G/3n*  -  G34>/3n')dS'  «  0.  (3.10) 

°  Jv  '  Jr 


Similarly,  it  can  be  shown  that 


A  +  u  G(iwP  -  Vx 

“  0  Jv  ~ 


:M) dV '  +  [  [A3G/3n*  -  G(n' *7’ )A]dV'  *  0. 
J  S 


(3.11) 


Next,  we  consider  an  actual  scattering  situation.  We  split  the  electro¬ 
magnetic  field  into  two  parts: 


*  -  AinC  .  *SC 

4>  *  4>  +  $ 


A  -  Ainc  +  ASC 


(3.12) 


SC  sc 

where  4>  and  A  satisfy  the  radiation  condition  at  infinity  and  the  differential 
equations 


n2.sc  ,  .  2. sc 
7  $  +  k  $  *0 


n2.sc  .  i  2. sc  n  .  |. 

V  A  +k  A  *0  in  Vj, 


_2  sc  .  .2 .sc  -1_  _ 

V  $  +  k  <t  *  e  V’P 

o  — 


V2ASC  +  k2ASC  -  P0(iu>P  -  VxM)  in  V. 


The  incident  field  satisfies  the  differential  equations 


(3.13) 


(3.14) 


_2.inc  .  2,inc  n 
V  4>  +  k  <I>  ■  0 


V2AinC  +  k2AinC  -  0  in  V  U  V, 


(3.15) 


We  now  go  on  to  apply  the  integral  formulas  (3.10)  and  (3.11)  to  $ 


and  A  .  Equation  (3.10)  gives 


4>SC  +  e"1  [  GV-PdV'  +  [  [($®C  -  4>SC)8G/3n' 

°  Jv  ~  J c  + 


-  G(3«®C/3n'  -  at^/Sn’ildS'  -  0,  r  G  V  (3 


where 


i>  (r)  ■  lim  *(r) 
r-*-S 


$_(r)  *  lim  4>(r), 
r-*-S 


On  the  surface  S  we  impose  the  boundary  conditions  that  *,  A  and  D’n  are 

xnc 

continuous.  Since  4>  is  continuous  everywhere  these  boundary  conditions 
imply  that 


-  <t!c 


3*?C/3n  -  3$SC/3n  -  c_1n-(P  -  P  )  -  e_1n-P 

+  -  o  —  — +■  - -  o - + 

where  we  have  made  use  of  the  fact  that  £  =  0  in  V^.  Suppressing  the  index  + 
we  arrive  at  the  following  expression 


*  Ainc 
9  *  4>  -  e 


-1  GV-PdV'  +  e"1  Gn’.PdS' 

1  Jv  "  °  is  ~  " 


and  equation  (3.17)  is  valid  at  points  inside  V. 
In  the  same  way,  equation  (3.11)  gives 


A  +  u  (iwP  -  V*M)dV'  + 

°  >V  ~  iS 


M)dV'  +  f  [(ASC  -  ASC 

is  “+ 


)3G/3n' 


-  G(n'.7)(A8c  -  ASC)]dS'  -  0,  r£V,  (3 


We  also  impose  the  boundary  conditions  that  A,  n-I5  and  n*H  are  continuous  on  S. 
Making  use  of  these  boundary  conditons,  equation  (3.5),  and  the  vector  formula 

nx(VxA)  *=  VA^  -  (n*V)A 

we  get 

A  =  AinC  +  wo  |  G(VxM  -  imP)dV'  -  uQ  J  G(n'xM)dS'.  (3.19) 

Equation  (3.19)  is  valid  at  points  inside  V.  In  the  surface  integral  of  equation 
(3.19)  we  have  used  the  definition 


M(r)  »  lira  M(r) . 
£-*-S 
reV 


To  sum  up,  we  have  the  following  set  of  equations  for  the  electromagnetic 
potentials,  valid  in  the  region  V, 


A  *  AlnC  +  p  f  G(VxM  -  iwP)dV'  -  p 

-  -  °  Jv  -  °  J 


G(n’xM)dS' 


.  .inc  -1 
$  *  $  -  c 

o 


GV-PdV'  +  e'1 
,  ~  o 


Gn'-PdS'. 


(3.20) 


Combining  equations  (3.2)  and  (3.20)  we  get  the  following  set  of  equations  for 
the  electromagnetic  field  quantities  E,  B,  a1 so  valid  in  V, 


_  _inc  ,  -ln 
E  *  E  +  e  V 
—  —  o 


GV-PdV'  -  E' 


-1V  I  Gn'-PdS' 
»  is - 


L 


+  iupo  |  G(VxM  -  iwP)dV'  -  iwpo  J  G(n'xM)dS’ 


B  -  Binc  +  p  Vx  G(VxM  -  iwP)dV'  -  p  Vx  G(n'xM)dS' 

--  O  Jv  ~  °Jc~- 


(3.21) 


where 


_inc  0.inc  .  .  ,inc 
E  ■  -  V$  +  iwA 
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and 


Binc  -  VxAinC. 


Employing  the  integral  formulas 


and 


to  equation  (3.21)  we  get 


GV-PdV'  - 

Jv  ~ 

Gn'-PdS'  +  V- 

s  J 

* 

GVxMdV’  * 

Jy 

Gn’xMdS'  +  7x 

s  J 

E  -  EinC  +  t 


^VV-f  GPdV’  +  u2u  [  GPdV'  +  iu>p  7x[  GHdV' 
’  Jv  ~  °  Jv  °  Jy 


B  •  BlnC  +  u  7x7x 
—  —  o 


[  GMdV*  -  luu  Vx[  GPdV ' , 

Jv  0  >V  “ 


Making  use  of  equation  (3.3)  and  the  facts  that 


Dinc  _  E  Einc 


(3.22) 


(3.23) 


(3.24) 


(3.25) 


and 


(7 


+  k2)  [  GPdV’  -  -  P 
JV  ~ 


(3.26) 


we  can  replace  the  set  of  equations  (3.24)  by  the  following  set  of  equations 


D  »  DinC  +  VxVx 


I  GPdV’  +  iioy  e  Vxf  GMdV' 

Jv  “  00  IV“ 

f  GPdV'. 
Jv 


B  -  BinC  +  u  VxVxf  GMdV'  -  iu>p  7x 
-  -  o  jv  - 


(3.27) 


From  now  on  we  restrict  our  calculations  to  linear  media,  i.e.,  there 
exist  £(r,k)  and  y(£,k)  such  that 
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(3.28) 


and  (3.32)  to  the  set  of  equations  (3.30)  we  arrive  at  the  following  set  of 
coupled  integral  equations  for  _P  and  M 


(g  +  2)/[3(g  -  1)]P  =  DlnC  +  (  7x7x(GP)dV’  +  ike-1  f  7x(GM)dV* 

(Jj  +  2)/[3(i3  -  1)]M  =  HinC  -  ike"1  [  7x(GP)dV’  +  f  7x7x(GM)dV' 

JV  }V 


(3.33) 


The  integral 


1  Vx7x(G^)dV' 

Jv 


(3.34) 


is  a  well-defined  quantity  for  any  point  inside  V  provided  that  ^  is  a  Htflder 
continuous  function  in  V.  This  can  be  seen  from  the  fact  that 


7xV*(G$)  =  [R  3£,(0'  ,$*)  +  S(R,e,^,)>2(£’) 


(3.35) 


as  R  »  |r_  -  _r'  i  0.  Here  0’  and  are  the  polar  angles  of  the  vector  R  *  £  -  jr' . 
Moreover,  F  and  |l  are  3x3  matrices  such  that 


|£(R,6\(J>,)|  -  0(R  )  as  R  0 


(3.36) 


and 

F(e',4>')sin  e ’ de Td<p*  -  0  (3.37) 

>S 

where  S  is  the  surface  of  the  unit  sphere.  The  principal  value  integral  (3.34) 
has  been  discussed  previously  in  the  literature  and  we  refer  the  interested 
reader  to  references  C5]  through  [83. 

The  singularity  inherent  in  the  set  of  integral  equations  (3.33)  makes 
it  difficult  to  determine  the  analytical  properties  in  the  complex  frequency 
plane  of  and  M  from  the  set  of  integral  equations  (3.33).  In  the  next  section 
we  will  start  with  the  set  of  differential  integral  equations  (3.27)  and  derive 
a  set  of  volume-surface  integral  equations  for  E  and  H.  This  set  of  integral 
equations  can  be  treated  by  means  of  the  Fredholm  determinant  theory. 
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IV.  Analytical  Properties  of  the  Field  Scattered  From  an 
Imperfectly  Conducting,  Finite  Body 

Based  on  the  differential  integral  equation  (3.27)  we  will  in  this  section 
first  derive  a  volume-surface  integral  equation  describing  scattering  from  an 
inhomogeneous,  imperfectly  conducting,  finite  body.  We  will  then  discuss  the 
analytical  properties  in  the  complex  frequency  plane  of  the  solution  of  this 
integral  equation. 

Employing  the  integral  formulas  (3.22)  and  (3.26)  to  the  set  of  equations 
(3.24)  we  arrive  at  the  following  set  of  volume-surface  integral  equations 


,inc 


Ck2(e  - 


1)GE  + 


VGg-1Ve ■ 


E]dV' 


(e  -  l)VG(n' -E)dS' 

Js 


+  ikZ 


(C  -  l)VG*HdV' 


(4.1) 


H  =  HlnC  + 


[k2(ti  -  1)GH  +  VGif ^iMGdV'  - 


(y  -  l)VG(n’ -H)dS’ 


-  ikZ-1  f  (e  -  DVGxEdV'. 

°  Jv 

In  deriving  the  set  of  equations  (4.1)  we  have  made  use  of  the  fact  that 


P  =  CQ(g  -  1)E  =  e_1 (c  -  1 )D, 


M  =  (p  -  1)H  -  U^lg"1(q  _  1}-' 


For  the  special  case  where  u  ■  1  and  o  =  0  the  set  of  equations  (4.1)  reduces  to 
the  integral  equations 


E  =  Einc  + 


[k2(g  -  1)GE  +  VGE-1Ve*E]dV'  - 


(e  -  l)VG(n' -E)dS' 


(4.2) 


H  -  HinC  -  ikZ-1 

-  -  o  j 


(c  -  l)VGxEdV 
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1 


[13] 

which  have  been  discussed  previously  in  the  literature 

Instead  of  considering  the  set  of  integral  equations  (4.1)  we  will  consider 
the  following  more  general  set  of  integral  equations 


E  -  [  Xt.-EdV  -  L.FdS'  -  f  M.-HdV'  -  ElnC, 
h  1  JS  1  h 

F  -  j  LjFdS '  -  |  L2-EdV'  -  j  M^HdV’  -  FinC, 

L_ •  HdV '  -  L  JdS'  -  f  M_-EdV’  *  Hinc, 
V  L  JS  i  JV 


H  - 


J  - 


L2JdS'  - 


r  6  V 


res 


r€V 


(4.3) 


L, • HdV '  -  f  M2 • EdV *  -  JinC,  r  6  S 

J  v 


where 


*  k 2 (l  -  1)GJ,+  g_1VGVe 
L_  -  k2(p  -  1)GI  +  a_1VGV{j 
Lx  '  (1  -  S)VG 

L2  -  k2(&  -  l)Gn  +  2_1 (n- VG)V| 
L3  -  (1  -  S)VG 

•  k2(p  -  l)Gii  +  y  *(n*VG)Vy 
Lj  -  (1  -  e) (n*VG) 

L2  -  (1  -  u)(n-VG) 

M.  •  -  ikZ  (D  -  1)VG* 

1  o 

«  -  ikz"1^  -  l)VGx 
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M,  =  ikZ  (Ci  -  1)  (n*VG) 
— 1  o  — 


Substituting 


M_  °  -  ikZ-1 (e  -  l)(n*VG), 
— z  o  — 


„inc  inc 

F  =  n*E  , 


(A. 4) 


Tinc  „inc 
J  =  n-H 


into  the  set  of  equations  (4.3)  one  can  verify  that  the  solution  of  the  set  of 
integral  equations  (4.1)  also  satisfies  (4.3).  In  fact,  having  the  solution  of 
(4.1)  we  only  have  to  substitute 


F  =  n*E, 

J  -  n-H 


(4.5) 


into  (4.3)  in  order  to  show  that  the  solution  of  (4.1)  satisfies  (4.3). 

Next,  we  will  show  that  the  solution  of  the  set  of  equations  (4.3) 

satisfies  the  set  of  equations  (4.1).  The  uniqueness  of  the  solution  of  the 

integral  equation  (4.2)  for  real  values  of  k  has  been  shown  in  reference  Cl3j. 

Without  going  into  the  details  of  the  proof  we  claim  here  that  the  uniqueness 

of  the  solutions  of  the  set  of  equations  (4.1)  and  (4.3)  can  be  shown  by  using 

the  technique  employed  in  reference  tl  1 3D •  Let  E^,  be  a  solution  of  (4.1) 

with  the  right  hand  sides  equal  to  E  and  H  .  Moreover,  let  E„,  F„,  H„,  J. 

— o  — o  — z  z  — Z  Z 

be  a  solution  of  (4.3)  with  the  right  hand  sides  equal  to  E  ,  F  =  n-E  ,  H  , 

— o  o  —  — o  — O 

Jq  =  n-]^.  Introduce  the  quantities 


E’  «  Ej  -  E2 


F '  =  n.EL  ~  F2 
H’  =  Hj  -  H2 
J'  -  n-Hj  -  J2 


(4.6) 
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Then  one  can  easily  see  that  E!' ,  F',  11 ' ,  J'  satisfy  the  set  of  homogeneous 
integral  equations 


E'  - 


F' 


H’  - 


L.-E'dV'  -  [  L.F'dS'  -  M-H’dV’  =  0 

v  1  's  1  h  1 

L  F'dS'  -  (  L  -E’dV’  -  [  M. *H'dV'  =  0 
s  1  J  v  i  iv  1 


(A. 7) 


L2*H'dV'  - 


L-J'dS '  -  M.-E'dV'  =  0 
S  Jv 


J'  -  |  L2J’dS'  -  |  JL^H'dV'  - 


M^E'dV'  *  0. 


From  the  uniqueness  theorem  it  follows,  however,  that  the  only  solution  of  the 
set  of  integral  equations  (4.7)  is  the  trivial  solution  E'  =  F'  *  H'  «  J'  =  0. 
Thus,  by  substituting  the  set  of  equations  (4.4)  into  the  set  of  integral 
equations  (4.3)  it  follows  that  the  solution  of  (4.3)  satisfies  (4.1)  for  real 
values  of  k.  From  the  principle  of  analytic  continuation  it  then  follows  that 
the  solution  of  equation  (4.3)  satisfies  equation  (4.1)  for  all  complex  values 
of  k. 

We  will  now  continue  with  the  solution  of  the  set  of  integral  equations 
(4.3).  First  we  will  transform  the  set  of  vector  integral  equations  (4.2)  into 
a  set  of  two  coupled  scalar  integral  equations.  For  that  reason  we  introduce 
a  set  of  volumes,  V^,  and  a  surface,  Sj  defined  by 


vl  *  t  :  £  *  +  r*.  *'€ 


{£.  :  L 


Lo 


r\  r'£  S}. 


Here  0  £  £  £  5  and 


r 

— o 


t  max  Ct'r"  +  dD 
r"«S 


where  t  is  an  arbitrary  unit  vector  and  d  >  0.  Notice  that  Vq  »  V  and  that  the 
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T 


regions  V  are  nonintersecting  regions  and  that  S  =  S  and  S,  are  two  non- 

o  1 

intersecting  surfaces.  Moreover,  define 


Ex(r), 


rev 


n(r)  = 


E/l-r^),  reVj 


Vi  “  2^),  r  G  V2 


3^),  r  G  V3 


«y(£-  iev4 


H,(r  -  5r  ),  r  e  Vc 


and 


a(r) 


F(r),  res 

—  o 

'  J(£  -  !<,>.  I^V 


The  set  of  integral  equations  (4.3)  can  then  be  transformed  into  the  following 
set  of  two  coupled  integral  equations 


n(r)  - 


^JnCr^dV’  - 


r  '  »y)o(£.'  )dS'  =  ninC(r);  r  G 


o(r)  -  j  r21(r,r,;Y)n(r,)dV’  -  j  ^(r.r '  ;Y)o(r  '  )dS’ 


(4.8) 


inc,  .  _  „ 

o  (r) ;  r£i 


Here 


5 

n  -  u  v  , 

i-0 


I  «  SoU  S1§ 
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f 


and  rn,  ri2,  r21,  r22  are  constructed  from  the  kernels  of  the  set  of  integral 
equations  (4.3)  in  the  same  way  as  the  quantity  T  used  in  section  IV  of  reference 
[2U.  Since  we  are  not  interested  in  the  actual  analytical  expressions  of  F^, 
ri2*  r21  and  ^22  we  leave  ic  to  the  interested  reader  to  find  these  expressions. 
From  the  definitions  of  the  kernels  of  the  set  of  integral  equations  (4.3)  we 
can  derive  the  following  asymptotic  expressions 

!rn<£  '  l')l  ~  |r  -  r '  -  (i  -  jVrJ"2 

when  r  £  V  and  r'  £  V.  and  r'  -*•  r  -  (i  -  1)r  ; 

—  i  —  j  —  —  -o 

Iri2(l  ~  r*)|  ~  |r  -  r'  -  (i  -  j)rj“2 

when  _r  £  and  _r '  £  and  _r’  -*•  _r  -  (i  -  j)^; 

lr21(I  ~  r')|  |r  -  r*  -  (i  -  j)rj'2 

when  r  £  S  and  r'  £  V.  and  r'  -*■  r  -  (i  -  j)r  ;  and 

—  l  —  j  —  —  — o 

1  r22(—  "  £.')  I  ~  \L  ~  L  ~  (i  -  jXrJ"1 

when  r  6  S,  and  r'  £  S.  and  r'  -*•  r  -  (i  -  j)r  . 

i  j  o 

Introducing  £  *  (r\,o)  the  set  of  integral  equations  (4.8)  can  be  cast 
into  the  following  operator  form 

(I  -  £)•!-  ilnc-  (4.9) 

The  operator  £  is  determined  by  the  kernels  F^,  r12>  r21,  r22  and  the  domains 

of  Integration  0,£.  It  is  easy  to  show  that  the  kernels  defining  the  iterated 
3 

operator  £  have  logarithmic  singularities  at  those  values,  of  (r,r')  where 
ri2*  r21  and  r22  are  sin8ular*  Since  any  solution  of  the  integral  equation 
(4.9)  also  satisfies  the  integral  equation 

<1  -  £3) •!  *  ijnc  (4.10) 
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where 


i}nc  =  (I  +  L  +  L2)*±inc 


we  will  consider  the  integral  equation 


(I_  -  £.  )  •§_  = 


(4.11) 


A  solution  of  the  integral  equation  (4.11)  is  derived  in  Appendix  C. 
Provided  that  e,  Ve/e,  p  and  Vjj/jj  are  analytic  functions  in  the  entire  complex 
frequency  plane  it  follows  from  the  analysis  presented  in  Appendix  C  that  poles 
are  the  only  singularities  of  the  operator  inverse  to  the  integral  operator 
defining  the  left  hand  side  of  (4.11).  The  location  of  these  poles  are  given 
by  those  values  of  y,  y  ,  where  there  exists  a  nontrivial  solution  of  the 

[2jn 

homogeneous  equation 


(i'£3>'i=0-  (4-12) 

To  conclude  this  section  we  have  shown  that  there  are  two  types  of 
singularities  in  the  complex  frequency  plane  of  the  field  scattered  from  an 
imperfectly  conducting,  inhomogeneous  body.  One  type  is  due  to  the  singularities 
of  the  incident  field.  The  other  type  is  due  to  the  scattering  body.  Sufficient 
conditions  for  the  latter  type  of  singularities  to  consist  of  only  poles  are  that 
the  scattering  body  is  of  finite  extent  and  that  e,  Vi/t,  u,  and  Vu/Ci  are  analytic 
functions  throughout  the  entire  complex  frequency  plane. 

The  approach  we  have  used  is  based  on  a  formulation  of  the  electromagnetic 

scattering  problem  in  terms  of  a  set  of  coupled  Fredholm  integral  equations  of 

the  second  kind.  This  approach  resembles  the  approach  used  in  deriving  the  so- 

[9-121 

called  Oseen  extinction  theorem  .  When  solving  for  the  natural  frequencies 

and  natural  modes  of  a  particular  scattering  body  it  might  be  advantageous  from 
the  numerical  point  of  view  to  use  other  formulations.  One  method  that  can  be 
used  is  an  integral  equation  of  the  first  kind  for  the  induced  polarization 
current  and  magnetization  current  used  in  reference  [17].  In  comparison  to 
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involving  perfectly  conducting  bodies,  the  formulation  used  in  this  note 
resembles  the  magnetic  field  formulation  whereas  the  formulation  used  in 
reference  £17]  resembles  the  electric  field  formulation.  Finially,  as  we 
have  seen  in  section  11,  scattering  from  a  homogeneous,  imperfectly  conducting 
body  can  be  formulated  in  terms  of  a  surface  integral  equation. 


3 


Y 


V.  Scattering  From  a  Perfectly  Conducting,  Finite  Body 
Within  a  Parallel  Plate  Region 

In  this  section  we  will  consider  electromagnetic  scattering  from  a  perfectly 
conducting,  finite  body  located  between  two  parallel  plates.  With  the  use  of 
image  theory  we  will  first  formulate  an  integral  equation  for  the  scattering 
problem  valid  for  real  frequencies.  Next,  we  will  find  the  analytic  continuation 
of  this  integral  equation  to  the  entire  complex  frequency  plane.  Finally,  we 
will  discuss  some  of  the  analytic  properties  of  the  solution  of  the  integral 
equation. 

The  geometry  of  the  scattering  problem  considered  is  depicted  in  figure  4. 
From  the  theory  of  images  it  follows  that  the  integral  equation  derived  from  the 
magnetic  field  formulation  can  be  cast  into  the  following  form: 

4  I  ~  "  I1"0-  (5.1) 

Here  is  the  induced  surface  current  density  on  the  scattering  object,  ^  is  the 
identity  operator,  L  is  an  integral  operator  defined  by 

L-i  =  f  n*0  (£*j.) DdS ' ,  (5.2) 

3S 

and  n  is  the  outward  unit  normal  of  the  surface  S  of  the  scattering  object. 
Moreover , 


£  -  £"  +  £+. 

G~  =  IG~, 

00 

G  (r,r' ;ik)  *  J  (4xRn>  1  exp(ikRn), 

•00 

r-  -  C(X  -  x’  -  2nd)2  +  (y  -  y')2  +  (z  -  z')2^, 


1^  *  A  ,  A  .  .  . 

_  -  -  xx  +  yy  +  zz, 


(5.3) 

(5.4) 

(5.5) 

(5.6) 

(5.7) 

(5.8) 
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(5.9) 


G+(r,r';ik)  =  Y  (4itR+)  ^  exp(ikR+), 

-  n  n 

—00 

Rn  *  +  x'  +  2nd)2  +  (y  -  y')2  +  (z  -  z')2]  ,  (5.10) 

and  d  is  the  spacing  between  the  perfectly  conducting,  parallel  plates. 

We  observe  that  the  series  representations  (5.5)  and  (5.9)  of  G  (_r,r] ;ik) 

and  G+(r,r';ik),  respectively,  converge  for  real  values  of  k  except  for  k  *=  mr/d 

[24] 

(n  is  an  integer)  .  We  will  later  return  to  the  question  of  convergence  of 
the  series  representation  of  the  Green's  function  in  connection  with  the  problem 
of  finding  representations  of  G"(r,r];v)  which  are  valid  in  the  entire  y-plane. 

We  will  now  investigate  the  kernel  of  the  integral  equation  (5.1), 

?*(£•!)  =  Vx(£~.j_)  +  Vx(G+.j_).  (5.11) 

Here  V  operates  on  the  first  argument  of  G“(r_,r_' ;ik) .  Making  use  of  Cartesian 
coordinates  (x,y,z)  one  has 

00 

?*(£*•!)  -  l  f(ir)2;-i  (5.12) 

—oo 

where 

f (R)  =  (4nR3)"1(ikR  -  l)exp(ikR) 


+ 

and  jT  have  the  dyadic  representations: 

D^Or.r  ’ )  =  -  (z  -  z’)xy  +  (y  -  y')xz  +  (z  -  z')yx  -  (x  -  x'  -  2nd)yz 

-  (y  -  y')zx  +  (x  -  x'  -  2nd)zy, 

D^(r,r')  =  -  (z  -  z')xy  +  (y  -  y*)xz  -  (z  -  z')yx  -  (x  +  x'  +  2nd)yz 

+  (y  -  y')zk  +  (x  +x'  +  2nd)2?. 

Thus,  L  can  be  split  into  two  parts  (c.f.  (5.2),  (5.3)  and  (5.12)), 

k  -  k  +  1?  (5.13) 

I 

i 
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+  t  m 

where  L”  are  two  integral  operators  with  kernels  K  ,  having  the  representations,  r 

00 

^(ryr'jik)  «  I  (£..£' Jik)  •  (5.14) 

—00 

Here 

K“(r,r*  ;ik)  -  f  (R*)r(r,r’), 

-  n 

and  are  two  dyadics, 

F~(r,r ' )  =  -  [(y  -  y')n  +  (z  -  z’)n  ]xx  +  (x  -  x'  -  2nd )n  xy 
“n -  y  z  y 

+  (x  -  x'  -  2nd)n  xz  +  (y  -  y')n  yx 
z  x 

-  C(z  -  z')n  +  (x  -  x'  -  2nd)n  3yy  +  (y  -  y')n  yz 

Z  X  Z 

+  (z  -  z' )n  zx  +  (z  -  z')n  zy 
x  y 

-  C(x  -  x'  -  2nd)n  +  (y  -  y')n  Uzz, 

x  y 

j£(r,£’)  -  L(y  -  y’)ny  +  (z  -  z’Jn^xx  +  (x  +  x’  +  2nd)nyxy 

+  (x  +  x'  +  2nd)n  xz  -  (y  -  y')n  yx 
z  x 

-  C(z  -  z')n  +  (x  +  x*  +  2nd)n  lyy  +  (y  -  y')n  xz 

Z  X  z 

-  (z  -  z')n  zx  +  (z  -  z')n  zy 

x  y 

-  C(x  +  x'  +  2nd)n  +  (y  -  y')n  Izz. 

x  y 

From  (5.14)  one  can  derive  the  following  asymptotic  expression  for 
|T(r,r’;ik), 

K^(r,r’;ik)  -  t  fan  1  exp(2ikd|n|)  +  Jgf (r, r')t>n  (5.15) 
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-2 

where  b  =  0(n  )  as  n  -*■  00  and  A  is  a  constant  dyadic.  Thus,  we  have 

n  - 

K'O^ir'jik)  =  ±  A  £  n  *  exp(2ikd|n|)  +  (r_,_r ' ;  ik)  (5.16) 

n^O 

where  R~  (1:,  r/ ;ik)  is  well  defined  for  all  real  values  of  k  and  £  +  r_’  .  Moreover, 
since 


\  n  *  exp(2ikdn)  «  -  ln[l  -  exp(2ikd)]  (5.17) 

n*=l 

+ 

the  series  representation  (5.14)  of  ](“ (r^, r_' ; ik)  converges  for  all  real  values 
of  k  such  that  k  #  mr/d.  Also  from  equation  (5.15)  the  series  (5.14)  converges 
conditionally  for  k  =  mr/d.  In  passing  we  want  to  point  out  that  the  series 

00 

K(r,r';ik)  =  \  [K^(r,r';ik)  +  K+(r,r ’ ;ik)]  (5.18) 

—  OO 

converges  absolutely  for  all  real  values  of  k.  This  can  be  seen  from  the  fact 
that  for  large  n  we  have  asymptotically: 

K“(r,r';ik)  +  K+(r,r';ik)  =  C(r,r')c  (5.19) 

=n -  -  - n 

-2 

and  c  =  0 (n  )  as  n 
n 

To  conclude  our  considerations  for  real  values  of  k  we  will  briefly 
investigate  the  convergence  of  the  series  representation  (5.4)  through  (5.10) 
of  the  dyadics  £~(jr,r/; ik) .  Employing  the  methods  used  when  deriving  (5.16) 
it  can  be  shown  that  one  representation  of  £“  (r  ,r_ ' ;  ik)  is  given  by 

00 

;ik)  =  £[( r,r/)  \  n"1  exp(2iknd)  +  £o(— *  ;ik>  (5.20) 

n=l 

+  + 
where  G~(_r,r/;ik)  is  well  defined  for  all  values  of  k  and  r_  ^  jr * ,  and  GjCr,ir') 

is  well  defined  for  _r  jt  jr' .  Thus,  from  equations  (5.17)  and  (5.20)  it  can  be 

seen  that  the  series  representations  (5.3)  through  (5.10)  of  £“(r,_r' ; ik)  diverge 

for  k  *  n-rt/d  (in  contrast  to  the  series  representation  of  JC“ (£_,£'  ;ik)) . 

Until  now  we  have  only  considered  the  operator  L  *  L(y)  when  y  is  a  purely 

imaginary  number.  In  order  to  find  an  analytic  continuation  of  £  from  the 
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imaginary  axis  into  the  complex  Y-plane  we  first  make  the  following  observation. 
The  series 


G^r.r’iY)  =  l  (AttR*)-1  exp(-yR*)  (5.21) 

a.  co 

converges  for  Re{y}  >  0.  Therefore,  the  definition  of  L(y)  from  G-(jr,jr';Y) 
by  equations  (5.2)  through  (5.10)  can  be  extended  to  all  complex  values  of  y 
such  that  Re{y)  i  0.  Next,  we  will  find  an  operator  ^(y) *  defined  almost 
everywhere  in  the  entire  complex  y-plane,  and  such  that  ^(y)  *  iXy)  for 
Re{y}  2:  0.  In  order  to  find  A  we  start  with  the  following  representation  of 


where 


L  *  L  + 

~~  — o 


l 

m=2,3 

£=0,1 


%  •!)(£> 


[^(r.r'jY)  +  KD(r,r,;Y)>f(r’)dS,, 


(5.22) 


r 


s  <Y;r,r' jm.OB^r.r/J-f  (r')dS', 


S"(y  *,r,r' ;m,£)  *  £  n^CR")"™  exp(-yR“). 

n=l  ‘ 

The  analytic  properties  in  the  complex  y-plane  of  the  functions  S”  (y;£,jt ' ;m,£) 
are  investigated  in  Appendix  A.  In  this  appendix  we  have  shown  the  existence 
and  given  a  method  of  constructing  two  functions,  E~(y;£,£' ;m,£),  defined  almost 
everywhere  in  the  entire  y-plane,  and  such  that 


E±(y;r,r,;m,l)  =  S* (y ;m,£)  (5.23) 

for  Re{y}  i  0.  Moreover  Z~(y£f£'  ;m,l)  are  analytic  functions  in  the  entire 
Y-plane  except  along  the  lines  Cn,  defined  by 

Cn  -  {Y  :  Re{y)  £  0,  Im{y)  =  mr/d,  n  integer}. 
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Thus,  we  can  define 


where  7.  (y ;  r  ,£'  ;m , '/.)  are  discontinuous  (see  figure  5). 
an  operator  =  My), 


A  =  L 
=  =o 


+  l 

m=2,3 

£=0,1 


(A'  +  A+  ) 

Tin  =m£ 


where 


S_(y  ;m,  (£,♦£.' )  *X(£.'  )dS ' 


(5.24) 


and  L  is  defined  by  equation  (5.22).  From  equations  (5.22)  through  (5.24)  it 
— o 

follows  that 


A(y)  =  L(y) 

for  Re(Y }  a  0.  Moreover,  A(y)  is  an  analytic  operator-valued  function  of  y 

except  on  the  lines  C  .  Of  course,  the  lines  C  can  be  considered  as  branch 

n  n 

cuts  in  the  complex  y-plane  of  the  operator  JL(y)  •  The  operator  A  is  the  analytic 
continuation  of  L  into  the  complex  Y-plane. 

We  now  go  on  to  consider  the  solution  of  the  integral  equation 

<y  I  "  A)  -i  “  IlnC  (5.25) 

for  any  complex  value  of  y.  Since  A(ik)  =  it  follows  from  the  uniqueness 

theorem  that  the  solution  of  the  integral  equations  (5.1)  and  (5.25)  are 
identical  for  y  *  ik.  Next,  we  will  consider  the  analytic  behavior  in  the  complex 

Y-plane  of  the  solution  of  the  integral  equation  (5.25).  It  is  easy  to  show 

2  2  f2l 

that  the  operator  A  »  defined  by  4  *f.  *  4*  (A'l.) »  is  °f  Hilbert  Schmidt  type  . 

The  integral  equation 

<{  I  -  42)-i  *  <j  I  +  4)*iinc  =  iinc 

r  2i 

can  then  be  solved  by  using  the  Fredholm  determinant  theory  .  Following  the 

procedure  in  reference  [2^  one  can  easily  see  that  the  inverse  operator 

(■|  I,  -  A)  *  has  two  types  of  singularities.  One  type  is  poles  at  those  values 
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of  y*  Yn»  where  there  exists  a  nontrivial  solution  of  the  homogeneous  integral 
equation 


C|  i  -  A<V>4  ■  K  *  °' 

The  other  type  of  singularities  of  the  inverse  operator  (4  ,1  -  4)  *  is  dis- 

1  2  " 

continuities  along  the  lines  C  ,  where  the  operator  t  I  -  A  is  discontinuous. 

n  2  1  -1 
For  the  case  where  all  the  poles  of  the  inverse  operator  (-  I  -  A) 

are  simple  poles  we  get  the  following  time-domain  representation  of  the  induced 

r  2~i 

current,  J(:r,  t),  due  to  a  delta-function  incident  TEM  wave. 


Hlnc  -  I  6(2  -  ct), 
—  — o 


I 

— o 


V 


J(r,t)  - 


U(ct  -  2o)[  l  <iinC<Yn)«p(Ynzo),hm>[<Bn.l^,]^)r1l^exp[Yn(c^  -  zq)1 

Ln,m 

00  fOD  -1 

+  (2vi)  1  l  +  iitn/d)exp(-ctc  +  iTrnct/d)]d£ 

-00  J  0  J 


Here 


,  inc 


±  (y)  =  nyy  exp(-Yz), 


-  dA/dy  evaluated  at  Y  =  y  , 

and  h^  are  linearly  independent  nontrivial  solutions  of  the  homogeneous 
integral  equations 


C|  i  -  A '  °- 

1  - 

where  hf  is  the  adjoint  operator  of  A_.  Moreover, 
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00  =  lim  ci  I  -  A(-5  +  in  +  iim/d)]-1 
^  n-Oi  2 

and  we  assume  that  the  inverse  operator  Cy  i  ~  4^)3  *  exists  as  y  approaches 

C  (from  above  or  below), 
n 

To  sum  up,  we  have  shown  that  the  singularities  in  the  complex  frequency 
plane  of  the  operator  inverse  to  the  integral  operator  of  the  magnetic  field 
formulation  are  poles  and  branch  cuts.  The  location  of  the  branch  cuts  is 
uniquely  determined  by  the  distance  between  the  parallel  plates.  Based  on  the 
theory  developed  here  we  might  in  the  future  undertake  a  numerical  study  of  the 
location  of  the  poles.  As  we  have  seen,  the  main  problem  is  to  find  a  represen¬ 
tation  of  the  integral  operator  that  is  valid  in  the  entire  complex  frequency 
plane.  There  still  remain  some  problems  to  be  solved  in  the  construction  of 
this  operator.  Needless  to  say,  the  "ordinary"  series  representation,  derived 
from  the  method  of  images,  is  not  valid  throughout  the  entire  complex  frequency 
plane. 

Based  on  eigenfunction  expansions  different  representation  of  the  Green's 
function  is  given  in  Appendix  B. 


Appendix  A 


In  this  appendix  we  will  consider  the  following  sums 


S(4;m,£)  =  l  n£Rnm  exp(-?Rn),  m  i  £ 


(A.  1) 


where 


2  L 

R  =  (n  +  an  +  0)  . 
n 

The  series  (A.l)  converges  for  Re{^}  >  0  when  m  and  £  are  any  finite  integers. 

It  also  converges  for  Re{?}  i  0  when  m  a  £  +  2.  In  this  appendix  we  will  find 
the  analytical  continuation  of  S(£;ra,£)  into  the  entire  C -plane  when  m  £. 

Let  N  be  a  finite  integer  such  that 

M  +  !  6 1  <  N  £  |a|  +  |b|  +  I.  (A. 2) 

We  then  have 


S(^;m,£)  =  H(?;m,£)  +  D(c;m,£) 


(A. 3) 


where 


and 


N-l 

H(c;m,£)=  £  n  R-"1  exp(-?R  ) 

n=l  n  n 


D(C;m, £)  =  Y  n^R  m  exp(-£R  ). 

■.1  n  n 

n=N 


Since  N  is  finite  it  easily  follows  that  H(s;m,£)  is  an  analytic  function  in  the 
entire  4-plane.  Next  we  consider  the  expansion 


Rnm  exp(-?Rn)  -  l  fk(4,m)n"‘m“k  exp(-nC).  (A. 4) 

k*0 


The  analytical  expression  for  f^U.m)  is  rather  complicated.  We  will  here  only 
consider  the  convergence  of  the  series  (A. 4).  We  have 
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1 


R_m 

n 


(A. 5) 


and  Che  expansion  (A. 5)  converges  for  n  i  N.  Moreover,  it  is  easy  to  show  that 

there  exists  a  finite  constant  A  such  that  R  >  nA  for  all  n.  We  then  have 

n 


|  R~m  exp(-cR  )|  <  (nA)  m|exp(-?n  £  a.  n”k) j 
n  n  k-0 

00 

*  (nA)  m|exp(-Cn  -  Ca,  -  A  l  a.n  k+1)| 

1  k=2  K 


*  (nA)  “ID  +  Cn(e)]exp(-Cn  -  tc^))  (A. 6) 


where 

C„(C> 

00 

-  1 

•  if "  -t+lv 

Furthermore,  we  have 

CO 

-  oo 

1  £  ® 

|cn(c>|*  l  klV:)-1 

«.=i 

I 

K|N"k+1(N/n)k_1  *  l  C|Nc|£(n:)_1n*fc  (A.7) 

_k=2 

k  J  £-1 

where 

c  - 

l  |a.|N-k+1  < 

Comparing  equations  (A. A)  and  (A. 7)  we  obtain 

|fk(C,m)|  <  CAm|Nc|k(k!)  1|exp(-Cn  -  Soijl,  k  2:  1.  (A. 8) 

From  equation  (A. 8)  it  follows  that  the  series  (A. 6)  converges  for  all  complex 
values  of  £.  In  passing  we  also  want  to  point  out  that,  from  the  construction 
of  fk(c,m),  it  follows  immediately  that  fk(e,m)  is  an  analytic  function  in  the 
entire  A-plane. 

From  equation  (A. 8)  it  also  follows  that 

00  oo  oo  00 

D(c;m,H)  =  l  l  f  (?,m)nl_,n  exp(-n;)  -  l  l  f.  (;,m)n?'“m  exp(-nc).  (A. 9) 

n-N  k-0  *  k-0  n-N 


I 


A2 


From  the  definition  of  the  r  function  it  follows  that  (c.f.  Appendix  D  of 
reference  C28]) 

tT*  exp (-nO  *  Zr(l)Tl  x*“  expC-(C  +  x)nDdx 

J  r\ 


(A. 10) 


and 


£  n  ^  exp(-ns)  *  Cr (£) D  *F  (?) 
n=N 


v-1. 


(A. 11) 


where 


MO  = 


x*  1  exp[-(x  +  OlG/Cl  -  exp(-£  -  x)3dx. 


(A. 11) 


Equation  (A. 11)  is  valid  for  Re  {5}  >  0  when  l  »  1  and  it  is  valid  for  Re{^}  £  0 
when  1^2.  Thus,  we  have 


D(c;m,*,)  *  l  fk(C,m)Fm_Jl(c)/r(m  *  *,) . 


k»0 


(A. 12) 


The  convergence  of  the  series  (A. 12)  follows  from  the  following  consideration: 
for  Rets)  >  0  there  exists  an  e  >  0  such  that 


ll  -  e-(C+%)l  >  c 


and  from  equation  (A. 11)  it  then  follows  that 

|F£(c)|  <  e*"1r(H)N_^|exp(-Nc) 


(A. 13) 


The  convergence  of  the  series  (A. 12)  is  now  clear  in  view  of  equations  (A. 8), 
(A. 11)  and  (A. 13). 

Introduce  the  function 


fl.1 

z  expC-(z  +  C)N]/[1  -  exp(-z  -  ?)3dz 
JC 


(A. 14) 


where  the  path  of  integration,  C,  is  along  the  real  axis  in  the  complex  z-plane. 
One  easily  sees  that  is  an  analytic  function  in  the  entire  c -plane  except 
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on  Che  lines  C 


f 


:  Re{?}  <  0,  Im{c)  ■  27rn,  n  integer}. 

We  also  notice  that 


*£(C  +  2irin)  -  *£(C) 

so  that  we  can  limit  our  investigation  of  $  k)  within  the  band 

{5  :  -  °°  <  Rek)  <  ®,  -it  £  Imk}  £  ir}.  From  equations  (A.  11)  and  (A.  14)  it 

follows  that 


*£k)  *  F£k)  for  Rek)  >  0. 

Thus,  the  function  <J^k)  is  an  analytic  continuation  of  the  function  F^( t). 
Making  use  of  complex  contour  integration  it  is  easy  to  show  that  the  limit 
values 


lim  *  (?  +  in),  5  <  0  (A. 15) 

n-+0± 

exist.  Moreover,  we  have 


lim  *  (?  +  in)  -  lim  +  in)  -  2vi(-0*  (A.  16) 

n-K)-  n-KH 

Thus,  we  can  define  a  function  Ak»m,£), 


Ak;m,!l)  -  l  f  k,»)*  k)/r(ra  -  £)  (A. 17) 

k-0  k  m~l 

and  for  Rek}  >  0  we  have 


A(c;m,£)  ■  D(c;m,£).  (A. 18) 

It  is  easy  to  show  that  the  series  (A.  17)  converges  for  c  <j£  Cfl  (c.f.  the  proof 
above  concerning  the  convergence  of  the  series  (A. 12)).  Thus,  the  function 
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A(c;m,£)  is  an  analytic  function  of  c  in  the  entire  ? -plane  except  on  the  lines 

C  .  Ue  use  A(£;m,£)  as  the  analytic  continuation  of  D(c;m,£)  into  the  entire 
n 

C-plane.  From  equations  (A. 3)  and  (A.  18)  it  then  follows  that  the  function 


I(C;m,£)  -  H(c;m,£)  +  A(?;m,£),  (A. 19) 

in  an  analytic  continuation  into  the  entire  c -plane  of  the  function  S(t;m,£) 
defined  by  equation  (A.l). 
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Appendix  B 


Using  eigenfunction  expansions  we  will  in  this  Appendix  find  a  representa¬ 
tion  of  the  scalar  Green's  function  different  from  the  one  used  in  section  V. 

The  Green's  function  Gj(i:,ir';Y)  of  a  parallel  plate  region  satisfies  the 
differential  equation 


E 


_  +  I2  +  L 

3p2  *  *  p2 


6(4,  .  4> * ) 6 (2:  -  2*)  (B.l) 


for  0  <  z  <  d  and  the  boundary  conditions 

3Gj 

— —  »  0  when  z  =  0,d. 


To  find  G^  that  satisfies  equation  (B.l)  we  proceed  in  the  usual  way.  Let 


OO  GO 


l  l  ^nn^ro^n^^m^n^  *  ^  *cos(nirz/d)cos(nirz  Vd)expCim($-$')3» 

ln=0  m=-°°  , 

P  >  P 

(B.2) 


*Y)  - 


OO  QO 


l  l  AnmIm^np^Im(^  P*)3  lcos(mrz/d)cos(nnz'/d)expCim ($-$’) 3, 
i  n=0  m«-<*>  , 

P  <  P 


where 

A*  ■  y2  +  (nir/d)2 

and  1^(0 »  Km(0  are  modified  Bessel  functions.  Note  that  (B.2)  satisfies  the 

boundary  conditions  (B.l).  To  determine  the  constants  A  ,  multiply  the 

run 

differential  equation  (B.l)  with  cos(rorz/d)  and  perform  the  integration 

rd  r2n  rp  '+E 

lim  (* • • *)pdpd$dz. 

e-*0  ao  ^0  ■'p'-E 

After  some  straightforward  algebraic  manipulations  combined  with  the  Wronskian 
for  the  modified  Bessel  functions  we  obtain 

Anm  ■  <Vd)'\(V',Km(V,)  (,-3) 
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where  e  ■  1  if  n  4  1  and  e  *  2.  Substituting  (B.3)  into  (B.2)  we  have 
no 


;y)  m  l  l  (en ifd)-1Im(Xnp<)Kn(Xnp>)cos(nnz/d)cos(nirz7d)expCim($-<t> '  )3 

n*°  (B.4) 


where  p<  (p. )  denotes  the  smaller  (larger)  of  p  and  p'. 

Next,  we  consider  the  analytical  properties  of  G^  in  the  complex  y-plane. 
In  doing  so  we  first  observe  that  the  series  representation  (B.4)  converges  for 
all  finite  values  of  y  such  that  Xfl(y)  t  0  when  £  +  jr * .  The  convergence  is 
easy  to  show  from  the  asymptotic  expansions  of  the  Bessel  functions.  From  other 
representations  of  G^(£,£'  ;y)  one  has 

»T)  ~  (4ir[r  -  jr’  |)_1exp(-y|r  -  r'  |)  as  £-+•£’. 

We  also  note  that  the  function  X  (y), 

n 

X  (y)  *  v'y^  +  (tvn/d)^ 
n 

is  multiple-valued  but  can  be  made  single-valued  by  introducing  branch  cuts 

parallel  to  the  negative  real  axis  and  starting  at  iiirn/d.  These  branch  cuts 

have  been  denoted  by  C  in  section  V.  (See  figure  5).  Since  X  =  0  at  iirn/d 

n  n 

it  is  clear  that  the  branch  cuts  for  X  also  make  I  (X  p)  and  K  (X  p)  unique 

n  m  n  m  n 

functions  of  y.  Thus,  Gj  (£,£’  ;y)  is  an  analytic  function  in  the  entire  y-plane 
except  on  the  branch  cuts  CR. 
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Appendix  C 


In  this  appendix  we  will  find  a  solution  of  the  integral  equation  (4.11) 
of  section  IV, 


(I  -  £3)-£  -  jj,. 


(C.l) 


The  quantities  F.  ,  £  and  ij^  have  been  defined  in  section  IV. 

The  operator  equation  (C.l)  has  the  following  explicit  representation 

n(r)  -  j  Aji (£»£.' ;y)n(r' )dV*  -  J  A 12 (£»£.* JY^Cr^dS'  *  nj(£), 

r  G  ft 


o(r)  - 


A2i(£«£' 5Y)n(i,)dV'  “  |  A22(r,£*;Y)a(£,)dS’  - 


(C.2) 


r  £  i 


Here  A^,  Aj2,  A22  have  a  logarithmic  singularity  at  those  points  where 

I’ll •  ri2,  r21,  r22  are  singular. 

The  method  we  will  use  in  solving  the  set  of  integral  equations  (C.2)  is 
based  on  the  Fredholm  determinant  theory.  Introducing  the  resolvent,  Rj  (£,£'  ;y)  , 
which  satisfies  the  integral  equation 


RjCr^'  ;y)  -  |^A11(£,£";Y)R1(£,',£,;Y)dV”  -  Au(£,£';y) 
we  can  transform  the  first  equation  in  (C.2)  into  the  following  equation 


n(r)  *  nj(r)  +  j  A12(£,£f  ;y)o(£' )dS'  +  j  ;Y)h1(£,)dV' 

+  f  Mi.i^YJA.oCr^YWr'OdV’dS" 
h*z  L 


(C.3) 

[18] 


(C.4) 


Substituting  the  expression  (C.4)  into  the  second  equation  of  (C.2)  we  arrive 
at  the  following  integral  equation  for  o(£). 
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o(r)  -  |  A' (£,r' ;Y)o(£*)dS’  -  o' (r) , 


r  e  E. 


Here, 

A'(r,r’;Y) 
AjCr.r’^) 
A^r.r'  »Y) 


o’(r) 


+  Aj^Cr.r'jy)  +  A^r.r’jy) 

f  A21(r,r";Y)A12(£",r' ;Y)dV" 

■'fi 

A„  (r.r’"  ;y)R.  (r,M  ,r";y)H,  ,(r",r'  ;Y)dV"’dV" 
Jfixn  21  ~~  1  -  -  12 

o^Cr)  +  a2(0  +  o3(r) 


(C  .5) 


c?2  (£.)  - 


a2i  (£.»£.’  ;y)n1(r.,)dv' 


a  (r )  -  [  A21(r,r,;Y)R1(r,,r";r)n1(r")dV’dV”. 

J  Jfixfl  4i 


All  the  operations  we  have  performed  so  far  are  strictly  formal  and  we 

will  now  go  back  and  examine  the  validity  of  each  operation.  First,  we  notice 

that  the  kernel  A^Or.r/  ;y)  is  a  square  integrable  function  on  the  domain 

flxfl.  Thus,  the  integral  equation  (C.3)  can  be  solved  by  using  the  Fredholm 

determinant  theory.  Moreover,  in  the  following  we  assume  that  g,  VI/ g,  p  and 

Vp/p  are  analytic  functions  of  y  *  s/c.  We  will  also  assume  that  the  incident 

field  is  an  analytic  function  of  y  implying  that  and  o^  are  analytic  functions 

of  y.  It  then  follows  from  equations  (4.3),  (4.8),  (4.10)  of  section  IV  and 

(C.2)  that  the  kernels  A^,  Aj2,  A^  and  A22  are  analytic  functions  of  y 

throughout  the  entire  complex  Y~plane.  Following  the  procedure  used  in  reference 

C2]  it  is  easy  to  show  that  the  resolvent,  R.  (ir,r/  ;y)*  is  an  analytic  function 

(1)  1 

of  y  except  for  those  values  of  y»  Y^  ,  where  there  exists  a  nontrivial 
solution  of  the  homogeneous  integral  equation 


n(r) 


Ai  i  <£.£'; 

n  11  n 


o. 


(C.6) 
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The  resolvent  has  a  pole  at  y  “  Yn  • 

Let  A  be  a  domain  in  the  complex  y-plane  such  that  Rj(r,jr';y)  is  analytic 

in  A.  (We  will  discuss  later  what  happens  at  the  poles  of  (r^r' ;y) ) .  From 

the  integral  equation  (C.3)  it  follows  that  R^(jt,j:';y)  has  a  logarithmic 

singularity  at  those  values  of  (£.£,')  where  Ajj  (£»£.'  ;y)  has  a  logarithmic 

singularity.  From  the  Schwartz  inequality  it  then  follows  that  Aj  and  A2 

are  square  integrable  functionson  0*0.  The  integral  equation  (C.5)  can  then 

be  solved  by  using  the  Fredholm  determinant  theory.  Since  o'  is  an  analytic 

function  of  y  on  A  it  follows  that  the  only  singularities  in  A  of  the  solution 

of  the  integral  equation  (C.5)  are  poles  at  y  *  y'  where  there  exists  a  non- 

n 

trivial  solution  of  the  equation 

o(r)  -  f  A' (r,r' ;y,)o(r,)dS'  =  0.  (C.7) 

-  JL  n  - 


It  now  follows  that  the  operator  inverse  to  the  integral  operator  defining 
the  left  hand  side  of  equation  (C.5)  has  two  type  of  singularities.  One  type  is 
poles  at  those  values  of  y  *  y^  where  there  exist  a  nontrivial  solution  of 
equation  (C.7).  The  locations  of  the  other  type  of  singularities  (if  any) 
coincide  with  the  locations  of  the  poles  of  the  resolvent  R^  (jr,r/  ;y) . 

In  order  to  investigate  the  behavior  of  the  solution  of  equation  (C.5) 
around  the  poles  of  R^  (ir.r/  ;y)  we  will  go  back  to  the  set  of  integral  equations 
(C.2).  Introducing  the  resolvent  Rj  (£.»£.' *y)  satisfying  the  integral  equation 


R2<£»l' »Y)  -  J  A22(r,r";y)R2(r",r’;y)dS"  *  A22(r,r* ;y) 


(C.8) 


we  have 


o(r) 


Oj(r) 


a2i  (£»£.' ;  y)  (i.f  >dv  *  + 


J  *2^-*-  ;Y)Oi(r.,)dS' 


+  f  R2(r,r,;Y)A„1(r,,r,,;y)n(r")dS'dV".  (C.9) 

Jz*n 

Here  n(r)  satisfies  the  integral  equation 


n(r)  -  A"(r,r’;y)n(r,)dV'  -  n"(r),  r£tl  (C.10) 


r 
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where 


A"(r,r';Y) 

A3(r,r';Y) 


AH  Y)  +  A3(r,r';Y)  +  A4(r,r';Y) 

(  A12(£,£";Y)A21(_r",r, ;Y)dS" 

■’Z 


A4  (£»£'  5y) 
n"(r) 


f  A, 2(r,r"’ ;Y)R  (r"’ ,r";Y)A21 (r",r’ ;Y)dV'"dV" 
Jexi 

+  n2(0  +  n3(r) 


n2(0 

n3(£) 


Ai  ;y)°i  (i_')ds' 

J  A12(r,r';Y)R2(r',r";Y)a1(r’,)dS'dS". 
IXI 


First,  we  observe  that  the  only  singularities  of  the  resolvent  R9(£,£' ;Y)  are 

(2)  Z 

poles  at  those  values  of  Y,  Yr  ,  where  there  exists  a  nontrivial  solution  of 
the  homogeneous  integral  equation 


o(r)  - 


A22(r,£,;Y^2))o(£,)ds’ 


0. 


(C.ll) 


Assuming  that  Y^^  for  all  possible  combinations  of  m  and  n  it  follows 

immediately  by  comparing  the  solutions  of  equations  (C.4),  (C.5),  (C.9)  and  (C.10) 
that  the  only  singularities  in  the  complex  Y-plane  of  the  solution  of  (C.2) 
are  poles. 

For  the  special  case  where  there  exist  an  m  and  an  n  such  that  Y^^  * 
we  return  to  the  set  of  integral  equations  (C.2).  Dividing  the  volume  ft  into 
equal  subvolumes  and  dividing  the  surface  Z  into  N2  equal  subsurfaces  we  can 
approximate  (C.2)  by  the  following  set  of  equations 


N 

-  I 

3-1 


Vj6j 


li 


(C. 12) 


51 


nOLj),  1  i  j  i  Nj 
a(r^  ) ,  Nj  +  1  S  j  «:  N 

IMI/Mj.  1  J=  j  ^  Nj 

6j  = 

j|r||/N2,  Nj  +  lSjiH 

jjn||  is  the  volume  of  the  region  Si  and  |jl||  is  the  area  of  the  surface  I.  The 
elements  are  determined  by  the  kernels  Ajj,  ^12’  ^21’  ^22: 

(r^r.^),  L±£  n»  IjG  n 

(-.Ij).  1*6  1.  ^G. 

<£!»£))»  Li^  £*  £j  G  E* 

The  solution  of  (C.2)  can  be  obtained  from  (C.12)  in  the  proper  limit  as  ® 

and  N2  -*■  «.  This  way  of  obtaining  the  solution  of  (C.2)  follows  closely  the 

Fredholm  method  of  obtaining  the  solution  of  an  "ordinary"  second  kind  integral 

equation^- From  the  solution  obtained  in  this  way  it  is  now  easy  to  show 

that  the  solution  of  (C.2)  is  an  analytic  function  of  y  except  at  the  poles 

[21 

where  there  exists  a  nontrivial  solution  of  the  equation 

(I  -  £3)-±  -  0. 

This  completes  our  proof  that  the  only  singularities  are  poles  of  the 

3 

operator  inverse  to  the  integral  operator  i  ~  • 


Here  N  *  +  N2> 
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